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Abstract 

We present results of a lattice computation of the matrix elements of the vector 
and axial-vector currents which are relevant for the semi-leptonic decays D — >■ K and 
D — > K*. The computations are performed in the quenched approximation to lattice 
QCD on a 24 3 x 48 lattice at = 6.2, using an (9(a)-improved fermionic action. In the 
limit of zero lepton masses the semi-leptonic decays D — >■ K and D — > K* are described 
by four form factors: f£, V, A\ and A2, which are functions of q 2 , where is the four- 
momentum transferred in the process. Our results for these form factors at q 2 = are: 
/+(0) = 0.67 + _ 7 8 , V(0) = 1.01 + _H , Ai(0) = 0.70 +^ , A 2 (0) = 0.66 +J° , which 
are consistent with the most recent experimental world average values. We have also 
determined the q 2 dependence of the form factors, which we find to be reasonably well 
described by a simple pole-dominance model. Results for other form factors, including 
those relevant to the decays D — > tt and D — >■ p , are also given. 
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1 Introduction 



Semi-leptonic decays of the heavy-light mesons have attracted considerable interest, as they 
play a crucial role in the determination of the elements of the Cabibbo-Kobayashi-Maskawa 
(CKM) mixing matrix and in the understanding of weak decays. In recent years, a machinery 
has been developed for calculating weak matrix elements from lattice simulations (for review 
lectures presented at recent lattice conferences see refs. [l]-[6]). D decays provide a good test 
of the method, since the relevant CKM matrix elements (V cs and V c d) are well constrained 
in the Standard Model. In addition comparisons between D and B decays reveal the size of 
non-leading terms in the Heavy Quark Effective Theory (HQET). 

The study of the decays D — > Kl + vi and D — > K*l + vi (and similarly D — > 7r/ + z^, D — > pl + vi) 
is particularly simple. They proceed via the spectator process in which a charm quark decays 
into a light quark (s or d) by emitting a H^-boson, which materializes into a lepton pair 
(/ + ,z^), as shown in Figure 1. With only a single hadron in the final state, there are no 
interfering diagrams or final-state interactions to take into account, unlike the situation in 
non-leptonic decays. 





u 



Figure 1: Feynman diagram relevant in semi-leptonic D — >■ K,K* decays. 

The non-perturbative strong interaction effects are contained in the matrix elements (K*\J IX \D) 
and {K\J^\D), where = 17^(1 — 75)0 is the relevant quark weak current 3 . In this paper we 
present the results of a lattice calculation of these matrix elements using the improved quark 
action proposed by Sheikholeslami and Wohlert [7]. We determine the dependence of the 
form factors on the momentum transfer (q), and study the phenomenological implications of 

3 This discussion applies equally well to the D — > it, p cases by modifying the appropriate flavour quantum 
numbers. 
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our results. Previous lattice studies of these decays, obtained using the Wilson quark action 
can be found in refs. [8]-[15], and using the Sheikholeslami-Wohlert action in ref. [16]. 

The plan of this paper is the following. In section 2 we review the experimental situation and 
give the general formulae necessary for the calculation of the D — > Kl + ui and D — > K*l + vi 
decay rates. In section 3 we describe the details of our simulation, and the methods used to 
determine the matrix elements (and hence the form factors) from the correlation functions 
computed on the lattice. The results, as is always the case in lattice simulations, are obtained 
for unphysically large values of the masses of the u and d quarks and have to be extrapolated 
to the chiral limit. Details of this extrapolation are presented in section 4, and in section 5 
we discuss the relation between the lattice vector and axial currents used in this study, and 
the corresponding continuum currents. In section 6 we present a compendium of all our 
results. Finally, in section 7, we study the implications of our results, comparing them with 
the experimental measurements from refs. [17]-[27], summarised in [28]-[29], and with other 
theoretical predictions [8]-[16] and [30]-[37]. 



2 Phenomenology 

Using Lorentz, parity and time-reversal invariance, the matrix elements for the decays 
D — » K and D — > K* , can be parametrized (in Minkowski space) in terms of invariant 
form factors as follows [9, 11, 12, 30] : 



(K\(V-AUD) = [p D + PK - q mh f + (q 2 )+ q ™ h 2 m W ) (1) 



(K;\(V - A),\D) = e* r %, (2) 



2V(q 2 ) 



+ i ~Z~T~Z — (Pd + PK*)»qp ~ i — —2m K *q f t(p D + p K *)p (3) 



m D + m K * 

; (PD+PK*)»qp ~ i — 5 

m D + m K * q 

where q^ = (pp — Pk(k*)Y is the four-momentum transfer, e* 13 is the polarization vector of 
the K* and / + '°, V, A and A\^ are dimensionless form factors. A can be written as 

A(q 2 ) = A (q 2 ) - A,(q 2 ) (4) 

2 m D + m K * 2 m D - m K * 2 

Mq) = 2m K * Mq) ~ 2m K * M<i ] (5) 

with A (0) = A 3 (0) and / + (0) = /°(0). In the limit of zero lepton masses, the terms 
proportional to f° in eq. (1) and to A in eq. (3) do not contribute to the total amplitude 
and hence to the decay rates. 



2 



The physical meaning of the different form factors is clear in the helicity basis, in which each 
of the form factors corresponds to a transition amplitude with definite spin-parity quantum 
numbers in the frame of the center of mass of the lepton pair. Pole dominance models [30] 
then suggest the following behaviour with q 2 : 



2^ ^(0) , ,_ 2 , A)(0) 



V (l ) = i 21 2 ' A o(Q / ~ 2/ 2 ' 

q 2 /mf- 1 — q /ttIq 

q j 



where mjp denotes the mass of the sc meson with spin J and parity P. This simple picture 
certainly has limitations. The pole-dominated form factor would vary very rapidly with 
q 2 near the end point. Another limitation is that for f°,Ai,A 2 and A 3 , the + and 1 + 
resonances are, in most cases, not known or only poorly established. On the lattice we can, 
in principle, determine form factors as functions of q 2 . Therefore, assumptions such as pole 
dominance are not needed. Indeed, an important motivation for lattice computations is the 
opportunity to test such assumptions from first principles. 

The total decay rates are given by: 

r<2 It/ 12 r(mn— rau ) 2 n 
G 2 XV I 2 r(m D -m K *) 2 . 

no^^ Vl) = |^ WW* 

x(|//+tf)| 2 +|/n4 2 )| 2 + l#V)| 2 ) W 

where X(q 2 ) = [m 2 D + m 2 K K , — q 2 ) 2 — 4m 2 D m 2 K K , . H° comes from the contribution of the 
longitudinally polarized K* and is given by [34] 

H°(q 2 ) = — Z^_{(m 2 D - m 2 K . - q 2 )(m D + m^A^q 2 ) - Am »f^ - A 2 (q 2 )} (10) 
2m K *\Jq lv m D + m K * > 



where px* is the momentum of the K* in the D-meson rest frame. H ± correspond to the 
contribution of the transverse polarizations of the vector meson and are given by [34] 



H±(q 2 ) = -{(m D + m^Mq 2 ) T 2m ^' V(q 2 )) (11) 

We now briefly summarise the experimental results for semi-leptonic decays of D-mesons, 
basing our discussion on the review articles [28]-[29]. The largest and best measured semi- 
leptonic decay is D — > Klv\. There have been several experiments ([17]— [21]) which have 
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measured the branching ratios B(D° — > K~l + vi) and B(D + — > K G l + ui). From these exper- 
iments and the total D° and D + lifetimes, one can calculate the D° and D + semi-leptonic 
decay rates, which should coincide by isospin symmetry. The world average value of the 
semi-leptonic width quoted in [29] is T(D — > Kl + ui) = (7.1 ±0.6) x 10 10 s _1 . However, 
by looking only at the D° — > K~l + v\ channel and assuming isospin symmetry, a different 
average value, T(D ->■ Kl + ui) = (9.0 ± 0.5) x 10 10 s _1 , is given in [28] 4 . 

The shape of the form factor is a measure of the decreasing overlap of the D and K wave 
functions as E K increases. CLEO has measured this shape with the largest sample of D° — > 
K~l + vi decays. Due to the phase space, the differential decay rate peaks at low q 2 . The 
factor |/x(<7 2 )| 2 increases with g 2 , changing by a factor of about 2 over the kinematical 
range of the decay. A good fit to the data is obtained using eq. (7) with a pole mass mi- = 
(2.00±0.12±0.18) GeV [17], which is in good agreement with the value of 2.1 GeV expected 
from the closest resonance with the proper quantum numbers, the D*. The measured value of 
the pole mass mi- agrees with earlier experiments but with a smaller error. /#(0) is obtained 
from the total semi-leptonic width, integrated over q 2 , by assuming pole dominance. The 
average value quoted in [29] is /-(0) = 0.70 ±0.03, whereas the average value quoted in [28] 
is /+(0) = 0.77 ± 0.04. 

The Cabibbo-suppressed decay D — > ttIui has also been observed. Since the ratio |Kd/K«| is 
known, assuming unitarity of the CKM matrix, from the comparison of the decays D — > irlvi 
and D — > Klvi it is possible to determine the ratio /+(0) //^(0) . This ratio is predicted 
theoretically to lie in the broad range 0.7-1.4. Mark III [20] gives a result of /+(0) //^(0) = 
1.0 +° 3 ± 0.1. In a recent analysis, CLEO gets a value /+(0) //± (0) = 1.29 ± 0.21 ± 0.11, 
[27]. The errors in this ratio of form factors are still very large. In addition in [28], a value 
for the rate T(D -)> ttZ+z^) of (1.2 ± 0.3) x 10 10 s' 1 is quoted. 

There have been a number of measurements of B(D — > K*l + ui), with both D° and D + 
mesons ([17],[21]-[25]). The average value of the width from these measurements is T(D — > 
K*l + vi) = (5.1 ± 0.5) x lO 1 ^- 1 [28]. A slightly different average value, F(D K*l + vi) = 
(4.5 ± 0.5) x 10 10 s _1 , is given in ref. [29]. The experimental results for the form factors 
V, Ai and A 2 are summarised in Table 1 and the results for the ratio of the decay rates of 
the longitudinal (H° contribution in eq. (9)) and transverse (H ± contributions in eq. (9)) 
K* are presented in Table 2. The total rate is dominated by the A\ form factor and the 
ratios of form factors are determined by fitting the angular distributions. 

For the (Cabibbo-suppressed) decay D — > plvi there only exists an upper limit for the 
branching fraction at 90 % confidence level B(D — > plv{) < 0.37 [20]. 

4 In [28] the most recent measurement of the CLEO collaboration ([17]) is included, whereas it is omitted 
in ref. [29]. 
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Exp. 


E-687 [22] 


E-653 [26] 


E-691 [23] 


World Ave. [28] 


V/A x 
A 2 /A 1 


1.74 ±0.27 ±0.28 
0.78 ±0.18 ±0.10 


2.00 ±0.33 ±0.16 
0.82 ±0.22 ±0.11 


2.0 ±0.6 ±0.3 
0.0 ±0.5 ±0.2 


1.90 ±0.25 
0.74 ±0.15 



Exp. 


A l 


A 2 


V 


World Ave. [28] 


0.61 ±0.05 


0.45 ±0.09 


1.16 ±0.16 



Table 1: Form Factors at q 2 = for the D — > K*l + vi decay. 



Exp. 


E-687 [22] 


E-653 [26] 


E-691 [23] 


World Ave. [29] 


r L /r T 


1.20 ±0.13 ±0.13 


1.18 ±0.18 ±0.08 


1.8 + °!±0.3 

-0.4 


1.2 ±0.1 



Table 2: Ratio of the longitudinal and transverse partial widths for the D — > K*l + ui decay. 

3 Details of the Simulation 

We work in the quenched approximation on a 24 3 x 48 lattice at (5 = 6.2, which corresponds 
to an inverse lattice spacing or 1 = 2.73 ± 0.05 GeV, as determined from the string ten- 
sion [38]. Other physical quantities lead to slightly different values for the lattice spacing 
(a -1 = 2.7 ~ 3.0 GeV) [39]. This uncertainty in the determination of the scale should be 
reflected in our results for dimensionful quantities. Our calculation is performed on sixty 
SU(3) gauge field configurations [38]. The gauge configurations and quark propagators were 
produced on the 64-node i860 Meiko Computing Surface at the University of Edinburgh. 
The SU(3) gauge fields were generated using the Hybrid Over- Relaxed algorithm, defined 
in reference [38]. The gauge configurations are separated by 2400 sweeps, beginning at con- 
figuration 16800. The quark propagators were calculated using an 0(a)-improved action 
proposed by Sheikholeslami and Wohlert, which we refer to as the SW-action [7], 

S S F W = S™-t^J2 q(x)F,AxKMx), (12) 

X,fJ,,V 

where S^f is the standard Wilson lattice action, 

Sf = E WM*) " « E [^)(1 " l,)U,(x)q(x + fi) + q(x + p)(l + lft )Ul(x)q(x)] }. (13) 

is a lattice definition of the field strength tensor and k is the hopping parameter. Periodic 
boundary conditions were employed in the spatial direction and anti-periodic in the temporal 
direction. The "improvement" is particularly important here since we are studying the 
propagation of quarks whose bare masses are around one third of the inverse lattice spacing. 
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The interpolating operators and currents which we use in this study are of the form: 

q^x) (1 + — r D)T\1- — r Dj q 2 (x) (14) 

where F is one of the 16 Dirac matrices. The matrix elements of these operators computed 
using the SW-action have no discretisation errors of O(a); the leading discretisation errors 
are of 0(a s a) [40]. 

Our statistical errors are calculated according to the bootstrap procedure described in 
ref. [38], for which the quoted errors on all quantities correspond to 68% confidence lim- 
its of the distribution obtained from 1000 bootstrap samples. 

We have computed light quark propagators at three values of the quark mass corresponding to 
k = 0.14144, 0.14226 and 0.14262, using an over-relaxed minimal residual algorithm with red- 
black preconditioning and point sources and sinks. The masses of the light pseudoscalar and 
vector mesons which are needed for this study were obtained in ref. [39] and are summarised 
in Table 3. Results extrapolated to the chiral limit (found to correspond to a hopping 
parameter K cr n = 0.14315 +2 ) and to the mass of the strange quark (k s = 0.1419 are 
also tabulated. Our D-mesons consist of a heavy quark (with k c = 0.129, where the subscript 
c stands for "charm") and one of the light antiquarks. The heavy quark with k = 0.129 has 
a mass approximately equal to that of the charm quark [41]. We use spatially-extended 
interpolating operators for the D-mesons (we use gauge-invariant Jacobi smearing on the 
heavy-quark field, described in detail in ref. [41]), but local operators of the form in eq. (14) 
for the light mesons. 

Further details on the lattice calibration, fitting procedures, mass spectrum and extraction of 
matrix elements of local operators between the vacuum and meson states, e.g. (M PS | QjsQ |0) 
can be found in references [39] and [42]. In Tables 3 and 4 we show a summary of the mass 
spectrum found which we will use below. In the following we only present those details of 
the calculation which are specific to semi-leptonic decays and cannot be found in the above 
references. 

In order to determine the matrix elements in eqs. (1) and (2) we compute the three-point 
correlation functions: 

C(p D: q : t D ,tj) = ^ e ^'V^(J B (^,f)J^(t J ,y)4(0,0)) (15) 

C^(p D ,q,t D ,tj) = ^e ifo -V^(J D (^,f)J^(^,j7)J?.(0,0)) (16) 

x,y 

where Jr> is a spatially-extended interpolating field for the D meson [41] and is an 
(3(a)-improved lattice operator corresponding to the continuum weak currents sj^(l — 7 5 )c 
or (77^(1 — 75 )c (q = u or d). Jk and Jf { * are local interpolating operators which can 
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meson 


1 meson 


0.14144 





14144 


0.298 


+ 2 
-2 


0.395 + I 

— 6 


0.14144 





14226 


0.259 


+ 2 

-2 


0.370 + t 

— 5 


0.14144 





14262 


0.241 


+ 2 
-3 


0.360 + l 

— b 


0.14226 





14226 


0.214 


+ 2 
-3 


0.343 + _\ 


0.14226 





14262 


0.192 


+ 3 

-3 


0.331 +jj 


0.14262 





14262 


0.167 


+ 3 
-4 


0.319 


k s =0.1419 




=0.14315 


0.181 


+ 9 
-8 


0.326 +jj 


K cr it=0. 14315 


^crit 


=0.14315 







0.290 



Table 3: Light-light meson masses in lattice units. For the pseudoscalar channel we fit over 
the time range t = 14 — 22. For the vector channel we fit over the time range t = 13 — 23. 





meson 


1 meson 


0.14144 


0.716 + _\ 




k s = 0.1419 


0.701 + _l 


0.732 +4 A 

-4 


0.14226 


0.692 + _\ 




0.14262 


0.683 + _\ 




K cr it=0- 14315 


0.665 +3 


0.697 +5 A 

-4 



Table 4: Heavy-light meson masses, k c = 0.129, in lattice units. We fit over the time range 
t = 11 - 22. 
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annihilate the light-light pseudoscalar and vector mesons respectively. The lattice vector 
and axial currents (eq. (14) with T = 7^ or j^j 5 ) are related to the continuum ones by 
renormalisation constants Z v and Z A ; we will discuss the determination of these constants 
in section 5. To evaluate these correlators, we use the standard source method [43]. We 
choose to = 24 (in lattice units) and symmetrize the correlators about that point using 
Euclidean time reversal [44]. The position of the light meson source is fixed at the origin 
and we have varied the time position of the current in the interval tj = 7 — 16. In Euclidean 
space, provided the three points in the correlators of eqs. (15) and (16) are sufficiently 
separated in time (tj,tn — tj 3> 1), the ground state contribution dominates and one finds 
for tj <t D : 



C»(p D ,q,t D = 2A,tj) -+ Zd }^\ Z * (l f? + % -Eo(Po)to e [E D{ f> D )-E K{ p D+mj 
y ' 2E D (p D ) 2E K (p D + q) 

A(i2)(K,p D + q\V»\D,p D )* (17) 

C^(p D ,q,t D = 24,tj) MM MM + ^ e -E D (p D )t De lE D (p D) -E K ,(p D+ ^ tj 

1 ' 2E D (p D )2E K ,(p D + q) 

A(t,)A( a )Y,tr(PD + q)(K;,p D + q\[V-Ar\D,p D y (18) 

r 

A(ct = 0) = -1 , A(a = 1,2,3) = 1 (19) 

E D is the energy of the D-meson and its wave-function factor, Z D (\p D \) = (0\J D (0)\D,p D ) , 
is a function of the meson momentum because we use spatially-extended interpolating oper- 
ators. E K (E K *) is the energy of the light pseudoscalar (vector) meson and the wave-function 
factors Z K (\p\) = (0\J K (0)\K,p) and Z K .(\p\) ((0\J%.(0)\K* r ,p) = e?(p)Z K .(\p\)) do not 
depend on the momentum of the meson (p ) because we have used local densities. The 
factors A(/i) in eq. (17) and A(/i)A(tt) in eq. (18) come from relating D-meson matrix ele- 
ments, which we obtain directly from the three-point correlation functions defined in eqs. (15) 
and (16), to those of the D-meson which we are interested in. 

The matrix elements have been computed for two values of the momentum of the D-meson 
((12a/ir)pD = (0,0,0), (1,0,0)); and all values of the momentum transfer q for which 
(12a/7r)|<f| < 2. In order to limit the systematic errors (and also statistical ones), we will 
only present results for matrix elements for which both the initial- and final-state mesons 
have three-momenta less than or equal to it /12a. To improve statistics, we average over all 
equivalent momenta and the different correlators, C^ ,a or C M , which lead to the same matrix 
element. 

The wavefunction factors and energies are obtained from fits to two-point correlation func- 
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Z 2 (p = (0,0,0)) 


z 2 D (P= is(i.o.o)) 


0.14144 
0.14226 
0.14262 


14.5 + * 

- 4 

12.6 + J 
12.0 + I 

— 5 


io.6 !j 

9.0 + * 
8.4 



Table 5: Pseudoscalar heavy-light meson wavefunctions (in lattice units), k c = 0.129. Fitting 
ranges are the same as those in Table 4. 

tions. At large times, t, the Euclidean correlators G 5 and Gij behave as follows: 

X 

-+ Z "(\p\)^y- cosh (^(P)[t--]J (20) 

Gij(t : p) = £ e ^(f,t)^(0,0)) 

X 

- H + f|)^ ( ' j7 " £ 9S !cosh ( £i ' (!7)[i -IO (21) 

where T is the length of the lattice in the time direction, P 5 and Vi are the pseudoscalar 
density and vector current with the appropriate flavour quantum numbers, and E§ and 
Ey the energies of the mesons. For light mesons we use continuum dispersion relations, 
i.e., E 5yV (p) = ^m\ v + p 2 and impose Z 5 (|p|) = Z 5 (|0|) [45]. These relations are well 
satisfied for momentum 7r/12a which is the highest one we have considered. For the D 
meson, as mentioned above, the wave-function factors Zb{p) depend on the momentum and 
it is necessary to fit the corresponding two-point correlators to the asymptotic expressions 
of eq. (20) not only for p= (0,0,0), but also for p = 7r/12a(l, 0, 0). We have constrained 
the energy to be E D (p) = \Jm 2 D + p 2 and therefore have performed only a one parameter 
fit in order to find Z D (\p\). The masses we have used in our study of semi-leptonic decays 
of D-mesons appear in Tables 3 and 4, whereas the wave-function factors appear in Tables 5 
and 6. 

Having determined the Z's and energies, all the factors multiplying the required matrix 
elements on the right hand sides of eqs. (17) and (18) are known, allowing us to determine 
the different form factors which appear in the matrix elements. The results presented in 
section 6 were obtained by fitting the different (/i, a) correlators, for each combination of 
quark masses and each momentum channel, to their respective asymptotic forms (eqs. (17) 
and (18)) in the time interval tj = 11 — 13. We have performed correlated fits, but we only 
allow for correlations between different timeslices (tj = 11, 12, 13) of a given (//, a) correlator 
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Zf(p= (0,0,0)) 


Z2(p= (0,0,0)) 


0.14144 


0.14144 


0.0081 +3 


0.0025 ^ 


0.14144 


0.14226 


0.0067 _J 


0.0021 +2 


0.14144 


0.14262 


_ „ „ „ ^ -1-4 

0.0062 , 

- 4 


0.0019 + _\ 


0.14226 


0.14226 


0.0056 


0.0017 


0.14226 


0.14262 


0.0052 + ^ 

- 4 


0.0015 



Table 6: Light-light meson wave-functions (in lattice units). Fitting ranges are the same as 
those in Table 3. 

at the same quark mass and momentum. 



4 Chiral Extrapolation 

We are interested in deriving the form factors for physical values of the charm, strange and 
light quark masses, for a range of values of the momentum transfer q. We obtain these 
by extrapolation from the three-point correlation functions of eqs. (15) and (16) computed 
with a fixed charm quark mass (corresponding to k c = 0.129) and for three values of the 
light quark mass (corresponding to Ki = 0.14144, 0.14226 and 0.14262) and two values of 
the strange quark mass (corresponding to ki s = 0.14144 and 0.14226). The extrapolation to 
the physical values of the light and strange quark masses proceeds as follows: 

i) For each set of three-momenta of the initial and final state mesons, we determine each 
form factor for the six combinations of light and strange quark masses. The masses 
are extrapolated to their physical values using: 

m D (K c ,Ki) = a PS + -^-(— — ) (22) 

m v (Ki,K 2 ) = a v + b v (^— + ^- — ) (23) 

ml s ( KuK2 ) = 6ps (-L + -L - JL) . (24 ) 

This extrapolation, together with the continuum dispersion relations, also determines 
the value of q 2 corresponding to each set of three-momenta for physical quark masses 
(see Tables 3 and 4). 

ii) For each momentum channel we extrapolate the form factors to the physical limit, using 
the full covariance matrix, assuming the following dependence on the quark masses: 



10 



by \ i , Am D (K c ,Ki) Amu ght (Ki,Ki s ) ( Am %M (K|, K h ) \ 2 

F{Ki s ,K t ) = a + b r + c + d (25) 

1TlD\K>ci Kcrit) m light[K C rit, K s ) \^light{K C rit: K s j J 

where AraD(K c ,Kj) and Amn g ht(ni, ni s ) (light stands for light pseudoscalar and light 
vector mesons) are defined by 5 

Am D (K c ,Ki) = m D (K c , Ki) - m D (K c , K crit ) (26) 

Am PS (Ki,Ki 3 ) = m PS {K h Ki s ) - m PS {K criu K s ) (27) 

Am v (Ki,Ki s ) = m v (ni,Ki a ) - m v (K crit ,K s ) (28) 

In the decay into vector mesons, we have not kept the quadratic term ( Amit aht( K iM s ) \ 2 
(from Table 3 it can be seen that, unlike in the case of the light pseudoscalar meson, 
this term is always smaller than 5% and has a negligible effect on the extrapolation of 
the form factors to the physical limit) and we end up with only three free parameters 
(a, 6, c). Thus, in the CT — > CT case, we fit the form factors to the following dependence 
on the quark masses: 



i 



F(K l .,K l ) = a + p(^-- J -)+7(- + 1 - — V + 5 (- + -- — ) (29) 
and in the CT — > 1 _ case we have assumed the following dependence: 

F( Kl . , Kl ) = a' + p (- - —) + i ( — + - - —) (30) 

\Kl K crit ) \Ki s Ki 



crit , 



where F represents a generic form factor. Note that, in contrast to some analyses (e.g. 
[13, 16]), we do not assume flavour symmetry between the active and spectator light 
quarks. Thus for example, the form factors extrapolated to the strange and critical 
quark masses, F K and F p , are: 



F K = « + 7 +S( (31) 

F p = a' (32) 



5 The dependence assumed in eq. (25) is motivated by the results of a Taylor expansion of g 2 (/c;,K;J 
around q 2 (K crit , k 3 ). 



11 



5 Renormalisation Constants Zy and Za 



In this section we discuss the difficulties in determining the form factors of semi-leptonic 
D — > K, K* decays, due to the presence of discretisation errors. Of course these errors 
are substantially reduced by the use of the improved action ([7], [40]), nevertheless even in 
this case we believe that, for currently accessible values of the lattice spacing, they lead to 
uncertainties of the order of 10% in the form factors. We will now attempt to justify this 
statement. 

The lattice currents and used in this study are related to the physical ones (V^ and 
A^) by renormalisation constants Zy and Za'- 

Z v Vf = V, ; Z A A\ = A, (33) 

When both quarks are light, at (5 = 6.2, Z v and Za are known to be about 0.83 and 1.05 
respectively [46]. For light quark masses the discretisation errors are very small, but for the 
charmed quark at (3 = 6.2 this is no longer the case. In previous simulations, using Wilson 
fermions, these effects were modelled by using an effective (mass-dependent) value of Zy and 
Za, or by assuming that the "conserved" vector current, i.e. the lattice current which would 
be conserved if the quarks were degenerate, is free of discretisation errors, [45]. However we 
wish to stress that the discretisation errors of 0(a s ma) and 0(m 2 a 2 ) are in general different 
for matrix elements of currents with different Lorentz indices and between different states. 
Thus they cannot be absorbed into an effective Z v or Za for all cases. To see this, note 
that there are discretisation errors due to the mixing of the currents with higher dimensional 
operators, e.g. the vector current can mix with aq\D li q 2 or a 2 qij fl D tl D tl q 2 . The behaviour of 
matrix elements of these operators with the external states is in general different from that 
of the currents. We have carried out an extensive study of these effects for the heavy-heavy 
vector current Q^Q 6 . 

Defining Zy^ by 

ff C 2 (t x ;p) p» 

Zv ~ Ct(t y ,t x ;p) E {M) 

where 

C 2 (t x ;p) = Y,e^ s (J P (x)4(0)) (35) 

X 

Cg(ty,t x ;0) = Y, e * S (Mx)V(y)4(0)) (36) 

and Jp and Jp are the interpolating operators which can create or annihilate the heavy-light 
pseudoscalar meson P. For degenerate quarks with k = 0.129, and using correlation functions 



3 Pull details of this study can be found in [47]. 
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with p = and \x = 4, we find Zy* = 0.9177 We note that this value differs by about 



2 



10% from that for Zy determined using light quarks and is a measure of the residual size of 
the discretisation errors using the improved action. Another important question is whether 
the effects are multiplicative. To study this we have computed Zy* for \Pi\ = \Pf \ = -p^, and 
using the current V 4 in the correlation function C 3 we find Zy* = 0.925(1) whereas using 
the current V 1 we find Zy* = 0.99(6). Unfortunately the latter error is too large for us to 
be able to determine whether the discretisation errors are the same in matrix elements 7 of 
V 4 and V 1 . In any case, because of hypercubic group invariance, there are no discretisation 
errors of order O(a) which affect spatial components of the currents in a different way to 
temporal components at each value of q 2 ; the leading discretisation errors for which this 
happens are of order 0(a 2 ) for the SW action (arising e.g. from matrix elements of the 
operator a 2 qi^D^D^q 2 ). 

For simulations using the Wilson fermion action, the discretisation errors are much larger. 
Kronfeld and Mackenzie ([49]) have argued that much of this uncertainty can be absorbed 
into a multiplicative m-dependent correction factor to the heavy quark propagator. Consider 
the continuum free propagator at zero three-momentum: 

J d 3 xS(x, 0) = iipV mt (37) 

Using Wilson fermions the analogous expression is 

l + 7o 



Jd 3 xS w (x,0) = -^p(l + m a)-^ (38) 

e -mt e -ma ^ 



1 TO _ — rat —ma 



where m = log(l + mofl) /a and mo is the bare mass. An important correction proposed by 
Kronfeld and Mackenzie is the e~ ma factor in eq. (39). With the SW action, even at the tree 
level, there can be no 0(ma) term in the correction factor and we find 

fd 3 xS sw (x,0) = (l + ](l + m a-—^ ))\l + m a)-^ (40) 

J 2 V 4 1 + mod J 

Numerically we estimate that this correction factor differs from 1 by only about 1.5%, for our 
charmed quark (re = 0.129), whereas we have seen that Zy* in heavy-to-heavy transitions 
is about 10% larger than Z v . This suggests that this correction factor accounts for only 
a modest part of the discretisation errors. This conclusion is supported by the behaviour 
of Zy* with mass, for which we find that the quadratic term (i.e. the 0(m 2 a 2 ) term) is 
relatively small [47]. 

Our conclusions from the above observations are as follows: 



7 There is some mild evidence [47] that Zy* may grow faster with mass for /j = 1 than for /i = 4 in 
heavy-to-heavy transitions. 
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i) The use of the free quark propagator is not a useful guide to the discretisation errors 
when using the SW action. For the Wilson action it is possible that the corresponding 
factor e~ ma accounts for part of the errors. However, the remaining uncertainties are 
not understood, and are in any case formally at least as large as for the SW-action (i.e. 
there are 0(a s ma) terms, and there is even no proof that there are no 0(ma) terms 
above tree level 8 ). 

ii) Formally there is no reason to believe that the discretisation errors can be absorbed 
into universal effective renormalisation constants Zy^ and Z e /^ . Even if 0(m 2 a 2 ) cor- 
rections are neglected, in which case the discretisation errors are independent of the 
Lorentz index of the current, these discretisation errors could be different for differ- 
ent form factors and they could have a different q 2 dependence than the form factors 
themselves. However, if for a given form factor, both the form factor itself and the 
discretisation errors have the same q 2 dependence, e.g. if the pole dominance formula 
is a good approximation to both, then the corresponding effective renormalisation con- 
stant, Zy f or Z e J* , is independent of q 2 (up to corrections of 0(m 2 a 2 )). Therefore 
if 0(m 2 a 2 ) corrections are neglected and assuming that the meson pole dominance 
model describes well the q 2 dependence of the different form factors and their discreti- 
sation errors, only the fact that discretisation errors could, in general, be different for 
different form factors, prevents the absorption of all of them into universal effective 
renormalisation constants. The situation is not improved by the use of "conserved" 
currents on the lattice. Here the situation is more difficult than in the evaluation of 
the Isgur-Wise function, where one only needs to evaluate a single form factor and this 
factorization is still possible. 

Given this discussion what can be done? 

In this study we recognise that discretisation errors are of 0(10%), and in spite of the 
discussion above, we assume that they can be modelled by Zy f and Z e / f (at least part of 
the errors can be so absorbed). Specifically for the vector current we take 

Z e v jS = 0.88 + _\ (41) 

which represents an approximate average of Z v = 0.83 (obtained with light quarks) and 
Zy* = 0.92 for the heavy-heavy current with the mass of the heavy quark corresponding to 
k = 0.129. For the axial current, a non-perturbative determination of Za (when both quarks 
are light 9 ) using a method based on Ward Identities [48], gives a value of Za = 1.05(1), [46]. 
A one-loop calculation in perturbation theory for the SW action ([50]) when the "boosted" 

8 Note that in nth order perturbation theory, in general, terms appear which behave as a n ma log™ (a) « ma 
[40]. 

9 In this case we do not have a non perturbative determination of Z e /* for the heavy-heavy current. 
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coupling suggested in [51] is used, gives Za = 0.97. Unlike the case of the vector current, 
perturbative and non-perturbative determinations of Z A do not agree for the light-light 
current, thus for the axial current we have decided to take 

Z e A ff = 1.05 +g (42) 

which corresponds to the non-perturbative determination of [46] , but with an increased lower 
error in order to account for the perturbative value mentioned above. 

Below we will discuss briefly the dependence of the results on the values in eqs. (41) and (42). 

Over the next few years, as high statistics simulations are performed at different values of 
P, it will become possible to study the discretisation errors in detail. 

6 Results 

In this section we present our results for the form factors obtained at six combinations of 
momenta of the initial and final state mesons, which are (in units of it /12a, and using the 
notation ->■ p K ,K*)- a) (0,0,0) ->■ (0,0,0), b) (0,0,0) (1,0,0), c) (1,0,0) (1,0,0), 
d) (1, 0, 0) (0, 0, 0), e) (1, 0, 0) (-1, 0, 0), and f) (1, 0, 0) ->• (0, 1, 0). The momenta of 
the initial state D-meson are fixed to be (0, 0, 0) or (1, 0, 0), but we average over all equivalent 
momenta of the light meson, so that, for example, case f) is really the average of the four 
terms in which the final state meson has momentum n/12a in the positive or negative y or 
z directions. 

The results for the form factors, together with the corresponding values of q 2 , are presented 
in Tables 7, 8 and 9. In these tables we also present the form factors extrapolated to physical 
quark masses for the decays D — > K,K*, following the procedure described in section 4. 

From the results of Tables 7, 8 and 9, one can in principle check the pole dominance relations 
given in eqs. (6) and (7). This is true in practice for some of the form factors. However, 
we have a very poor determination of both, the scalar (0 + ) and axial-vector (1 + ) meson 
masses. Thus, we have decided to extract both the pole masses (rrijp) and the form factors 
at q 2 = 0, by fitting the chirally-extrapolated data to the pole dominance model. In the 
case of f + ,V and A we will compare the masses of the 0~ and 1~ mesons obtained from 
the pole dominance fit with the lattice results obtained in our simulation (Table 4). 

6.1 The exclusive 0~ — > 0~ case: 

With the method described in the former sections we have found for the D — > K decay, 
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V L) 








q 2 a 2 


f + (Q 2 )/Zy ff 


f°(Q 2 )/Zy ff 


(0,0,0) 


(0,0,0) 


0.14144 





14144 





175 


+ 2 

- 2 






1.03 


+ 4 

-4 











14226 





209 


+ 2 

-2 






0.99 


+ 5 

-5 






0.14226 





14144 





188 


+ 2 
-2 






1.04 


+ 5 
-5 











14226 





229 


+ 3 

-3 






1.00 


+ 6 

-6 






0.14262 





14144 





195 


+ 3 

-3 






1.05 


+ 6 

-6 











14226 





241 


+ 4 
-3 






1.03 


+ 8 
-7 










K s 





235 


+ 8 

-9 






1.01 


+ 8 

-7 


(0,0,0) 


(1,0,0) 


0.14144 





14144 





033 


+ 1 

- 1 


0.91 


+ 3 

-3 


0.87 


+ 3 

-3 











14226 





052 


+ 1 
- 1 


0.88 


+ 4 
-3 


0.81 


+ 3 
-3 






0.14226 





14144 





037 


+ 2 
-2 


0.90 


+ 4 
-4 


0.85 


+ 3 
-3 











14226 





057 


+ 2 

-2 


0.86 


+ 5 

-4 


0.78 


+ 4 

-3 






0.14262 





14144 





039 


+ 2 
-2 


0.89 


+ 4 
-4 


0.83 


+ 4 
-4 











14226 





060 


+ 3 
-2 


0.84 


+ 6 
-5 


0.75 


+ 5 
-4 






^crit 









052 


+ 4 

-4 


0.88 


+ 5 

-5 


0.80 


+ 4 

-4 


(1,0,0) 


(1,0,0) 


0.14144 





14144 





134 


+ 1 

- 1 


1.13 


+ 7 

-7 


0.91 


+ 6 

-7 











14226 





155 


+ 2 
- 1 


1.12 


+ 8 
-8 


0.82 


+ 7 
-7 






0.14226 





14144 





138 


+ 2 

-2 


1.09 


+ 11 

-11 


0.77 


+ 10 

-11 











14226 





162 


+ 2 

-2 


1.07 


+ 13 

-14 


0.64 


+ 11 
-11 






0.14262 





14144 





141 


+ 2 
-2 


1.00 


+ 14 
-15 


0.58 


+ 14 
-15 











14226 





165 


+ 3 
-2 


0.94 


+ 19 

-21 


0.40 


+ 16 

-17 






^crit 









157 


+ 4 

-4 


1.16 


+ 20 

-19 


0.73 


+ 15 

-15 



Table 7: Form factors for 0~ — > 0~ decay with momenta in units of ^ and k c = 0.129. For 
f°(q 2 ), the channel (1,0,0) — > (1,0,0) has not been considered in the pole dominance fit, 
because we feel we do not control its chiral extrapolation. 
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Pd 


Pk 






q 2 a 2 


f + (Q 2 )/Zy ff 


f°(q 2 )/Z e v ff 


(1,0,0) 


(0,0,0) 


0.14144 





14144 


0.147 


+ 2 

-2 


1.25 


+ 5 

-5 


0.98 


+ 4 

-4 











14226 


0.185 


+ 2 
-2 


1.33 


+ 7 
-7 


0.96 


+ 4 
-4 






0.14226 





14144 


0.163 


+ 3 
-2 


1.27 


+ 7 
-7 


0.99 


+ 5 
-4 











14226 


0.208 


+ 3 

-3 


1.40 


+ 10 

-10 


1.00 


+ 6 

-5 






0.14262 





14144 


0.172 


+ 3 
-3 


1.26 


+ 9 
-10 


1.00 


+ 5 
-5 











14226 


0.223 


+ 4 

-3 


1.43 


+ 14 

-15 


1.01 


+ 6 

-7 






^crit 






0.217 


+ 9 

- 9 


1.48 


+ 16 

-17 


1.06 


+ 7 

-8 


(1,0,0) 


(-1,0,0) 


0.14144 





14144 


- 0.140 


+ 1 

-1 


0.60 


+ 3 

-3 


0.71 


+ 4 

-3 











14226 


- 0.119 


+ 2 
-1 


0.58 


+ 3 
-3 


0.69 


+ 4 
-4 






0.14226 





14144 


- 0.136 


+ 2 

-2 


0.59 


+ 5 

-5 


0.71 


+ 6 

-5 











14226 


- 0.113 


+ 2 

-2 


0.58 


+ 6 

-5 


0.69 


+ 7 

-6 






0.14262 





14144 


- 0.133 


+ 2 
-2 


0.58 


+ 7 
-7 


0.69 


+ 9 
-8 











14226 


- 0.109 


+ 3 

-2 


0.57 


+ 10 

- 8 


0.67 


+ 11 

-10 






^crit 






- 0.117 


+ 4 
-4 


0.63 


+ 8 
-8 


0.73 


+ 9 
-9 


(1,0,0) 


(0,1,0) 


0.14144 





14144 


- 0.003 


+ 1 
-1 


0.78 


+ 3 
-3 


0.79 


+ 3 
-3 











14226 


0.018 


+ 2 
-1 


0.75 


+ 4 
-4 


0.73 


+ 4 
-4 






0.14226 





14144 


0.001 


+ 2 

-2 


0.76 


+ 4 

-4 


0.76 


+ 4 

-4 











14226 


0.025 


+ 2 
-2 


0.73 


+ 5 
-5 


0.70 


+ 5 
-5 






0.14262 





14144 


0.004 


+ 2 
-2 


0.74 


+ 6 
-6 


0.74 


+ 5 
-5 











14226 


0.028 


+ 3 

-2 


0.71 


+ 7 

-7 


0.67 


+ 7 

-7 






^crit 






0.020 


+ 4 
-4 


0.75 


+ 5 
-6 


0.72 


+ 5 
-5 



Table 7: (cont.) 
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Vn 

f L) 








q 2 a 2 


V(q 2 )/Z e v ff 


(0,0,0) 


(0,0,0) 


0.14144 






14144 
14226 


0.103 
0.119 


+ 2 
-2 
+ 3 
-3 










0.14226 






14144 
14226 


0.103 
0.121 


+ 3 

-3 
+ 5 

-5 










0.14262 






14144 
14226 


0.104 
0.124 


+ 4 
-4 
+ 7 

-6 










^crit 






0.112 


+ 5 

-5 






(0,0,0) 


(1,0,0) 


0.14144 





14144 


- 0.010 


+ 2 

- 1 


1.50 


+ 5 

- 5 











14226 


0.000 


+ 2 
-2 


1.47 


+ 7 
- 7 






0.14226 





14144 


- 0.012 


+ 2 

-2 


1.45 


+ 8 

- 9 











14226 


- 0.001 


+ 3 
-3 


1.41 


+ 10 
-11 






0.14262 





14144 


- 0.012 


+ 3 
-2 


1.43 


+ 12 
-12 











14226 


- 0.001 


+ 4 

-4 


1.36 


+ 16 

-15 






^crit 






- 0.009 


+ 3 
-3 


1.40 


+ 13 
-12 


(1,0,0) 


(1,0,0) 


0.14144 





14144 


0.083 


+ 2 
-2 


2.7 


+ 4 
- 4 











14226 


0.095 


+ 3 
-2 


2.9 


+ 5 
- 6 






0.14226 





14144 


0.082 


+ 3 

-2 


3.4 


+ 7 

- 7 











14226 


0.095 


+ 4 
-3 


4.1 


+ 9 
- 9 






0.14262 





14144 


0.082 


+ 3 

-3 


3.7 


+ 10 

-10 











14226 


0.096 


+ 5 

-5 


5.2 


+ 13 

-13 












0.086 


+ 3 
-4 


3.7 


+ 8 
- 7 



Table 8: Vector form factor (V(q 2 )) for 0~ — > 1 decay with momenta in units of ^ and 
k c = 0.129. The channel (1,0,0) — > (1,0,0) has not been considered in the pole dominance 
fit, because we feel we do not control its chiral extrapolation. 
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Pd 


Pk* 






q 2 a 2 


V(q 2 )/Z e v ff 


(1,0,0) 


(0,0,0) 


0.14144 





14144 





066 


+ 3 

-3 


i.46 +- 











14226 





085 


+ 4 
-4 


1.39 +_\\ 






0.14226 





14144 





068 


+ 4 
-3 


1.36 ™ 











14226 





088 


+ 6 

-5 


1.30 + _H 






0.14262 





14144 





069 


+ 5 
-4 


1A2 +28 











14226 





092 


+ 8 

-7 


1A2 +38 






^crit 









079 


+ 5 

-6 


1.09 + \l 

— 16 


(1,0,0) 


(-1,0,0) 


0.14144 





14144 


- 


191 


+ 2 

-2 


0.90 + _\ 











14226 


- 


179 


+ 3 
-2 


0.90 + I 






0.14226 





14144 


- 


192 


+ 3 

-2 


0.85 !» 











14226 


- 


179 


+ 4 

-3 


0.87 + )l 

— 15 






0.14262 





14144 


- 


192 


+ 3 
-3 


o.8i 











14226 


- 


179 


+ 5 

-5 


0.81 + _H 






^crit 






- 


188 


+ 3 
-4 


0.80 +J3 


(1,0,0) 


(0,1,0) 


0.14144 





14144 


- 


054 


+ 2 
-2 


1.22 +6 fi 

— b 











14226 


- 


042 


+ 3 
-2 


1.18 +1 






0.14226 





14144 


- 


055 


+ 3 

-2 


1.11 + 8 

-11 











14226 


- 


042 


+ 4 
-3 


1.06 






0.14262 





14144 


- 


055 


+ 3 
-3 


1.00 ™ 











14226 


- 


042 


+ 5 

-5 


0.90 






^crit 






- 


051 


+ 3 
-4 


1.02 + ] 2 A 

-14 



Table 8: (cont.) 
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Pd 
(0,0,0) 



Pk 
(0,0,0) 



0.14144 



0.14226 



0.14262 



0.14144 
0.14226 
0.14144 
0.14226 
0.14144 
0.14226 



0.74 
0.70 
0.76 
0.72 
0.78 
0.73 
0.76 



+ 2 
-2 
+ 3 
-2 
+ 3 
- 3 
+ 4 
-4 
+ 5 
-5 
+ 6 
-6 
+ 4 
-4 



Mq 2 )/z_ 



Mq 2 )/z_ 



eff 



(0,0,0) 



(1,0,0) 



0.14144 



0.14226 



0.14262 



^crit 



0.14144 
0.14226 
0.14144 
0.14226 
0.14144 
0.14226 



0.70 
0.65 
0.71 
0.66 
0.71 
0.64 
0.72 



+ 2 
-2 
+ 2 
-2 
+ 3 
-4 
+ 4 
-4 
+ 5 
-5 
+ 6 
-6 
+ 4 
-4 



0.59 
0.51 
0.63 
0.53 
0.61 
0.49 
0.69 



- 7 
+ 7 

- 7 
+ 9 

- 9 
+ 11 
-12 
+ 13 
-14 
+ 17 
-18 
+ 10 
-10 



0.73 
0.72 
0.73 
0.72 
0.74 
0.72 
0.72 



+ 2 
-2 
+ 3 
-3 
+ 3 
-3 
+ 4 
-4 
+ 5 
-5 
+ 6 
-5 
+ 4 
-4 



(1,0,0) 



(1,0,0) 



0.14144 



0.14226 



0.14262 



^crit 



0.14144 
0.14226 
0.14144 
0.14226 
0.14144 
0.14226 



0.76 
0.71 
0.66 
0.60 
0.49 
0.41 
0.68 



+ 7 
- 7 
+ 9 
-9 
+ 11 
-12 
+ 15 
-16 
+ 19 
-21 
+ 26 
-27 
+ 14 
-15 



1.4 +n 

-li 

1.5 ^ 
0.7 ^ 
0.7 + 20 

-20 

0.7 + 26 

-26 

9 +30 

U-y -30 

1.4 +23 

-23 



0.92 
0.96 
0.74 
0.78 
0.53 
0.57 
0.78 



+n 

-li 

+ 15 
-15 
+ 18 

-19 
+ 25 
-25 
+ 27 
-30 
+ 36 
-38 
+ 23 
-24 



Table 9: Axial form factors for 0~ — > 1~ decay with momenta in units of ^ and k c = 0.129. 
For A 2 , the channel (1,0,0) — > (1,0,0) has not been considered in the pole dominance fit, 
because we feel we do not control its chiral extrapolation. 
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Pd 


Vk 






Mi 2 )/z e A ff 






(1,0,0) 


(0,0,0) 


0.14144 





14144 


0.69 


+ 3 
-3 


0.67 


+ 39 

-38 


0.88 


+ 6 

- 7 











14226 


0.66 


+ 3 
-3 


0.93 


+ 50 
-46 


0.88 


+ 9 
- 9 






0.14226 





14144 


0.69 


+ 4 
-4 


0.57 


+ 55 
-52 


0.92 


+ 9 
-10 











14226 


0.66 


+ 5 

-4 


0.92 


+ 78 
-71 


0.91 


+ 13 

-14 






0.14262 





14144 


0.69 


+ 6 
-6 


0.37 


+ 75 
-71 


0.99 


+ 12 
-14 











14226 


0.65 


+ 7 

-6 


0.76 


+100 
-100 


0.99 


+ 18 

-21 






^crit 






0.66 


+ 5 

-5 


0.51 


+ 75 

-66 


0.90 


+ 12 

-14 


(1,0,0) 


(-1,0,0) 


0.14144 





14144 


0.57 


+ 4 

-4 


0.36 


+ 8 

- 7 


0.44 


+ 3 

-3 











14226 


0.53 


+ 5 
-5 


0.30 


+ 9 
- 8 


0.44 


+ 3 
-4 






0.14226 





14144 


0.59 


+ 7 

-8 


0.37 


+ 12 

-12 


0.45 


+ 4 

-5 











14226 


0.56 


+ 9 

-9 


0.31 


+ 15 

-15 


0.47 


+ 5 

-6 






0.14262 





14144 


0.58 


+ 12 
-12 


0.34 


+ 19 
-19 


0.46 


+ 7 
-8 











14226 


0.58 


+ 13 

-15 


0.30 


+ 23 

-23 


0.50 


+ 8 

-9 






^crit 






0.58 


+ 8 
-9 


0.36 


+ 15 
-15 


0.45 


+ 5 
-6 


(1,0,0) 


(0,1,0) 


0.14144 





14144 


0.61 


+ 2 
-2 


0.40 


+ 6 
- 6 


0.63 


+ 3 
-3 











14226 


0.55 


+ 3 
-3 


0.34 


+ 8 
- 7 


0.61 


+ 3 
-3 






0.14226 





14144 


0.59 


+ 4 

-4 


0.34 


+ 9 
-10 


0.63 


+ 4 

-4 











14226 


0.53 


+ 5 
-5 


0.28 


+ 12 
-12 


0.61 


+ 5 
-6 






0.14262 





14144 


0.54 


+ 6 
-7 


0.23 


+ 15 
-17 


0.64 


+ 7 
-6 











14226 


0.49 


+ 8 

-8 


0.20 


+ 20 

-20 


0.61 


+ 8 

-8 






^crit 






0.60 


+ 5 
-5 


0.36 


+ 12 
-12 


0.59 


+ 6 
-6 



Table 9: (cont.) 
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0.67 ^[a" 1 ] 

0.91 ^[a- 1 ] (43) 

The pole mass , agrees reasonably well with the value 0.732 ^ [a -1 ] quoted in Table 4 
and the experimental value (2.00 ± 0.12 ± 0.18) GeV [17] quoted in section 2. The result 
found for is also consistent with the value of 2.3 ± 0.2 GeV obtained 10 in the lattice 
simulation of ref. [13]. We also see /^(0) is equal to /^(0) within the errors. These values 
agree with, and update the preliminary results presented in [52]. In Figs. 2 and 3 we show 
the form factors / + and f° as a function of q 2 for two combinations of the light quark masses: 
Ki = Ki s = 0.14144 and Ki = K cr i t , Ki 3 = k s . In all cases the solid line corresponds to the 
comparison of the pairs (q 2 Jk'°(q 2 )) with the (^-dependence of the form factors determined 
from a two-parameter pole dominance fit to our data (giving the parameters in eq. (43), 
for the case ki = n crit , ki s = k s ). Crosses correspond to the form factors at q 2 = (up to 
a factor Zy*) determined in this way. In the case of / + , we also compare (dashed lines) 
the (^-dependence of our data with that determined from a one-parameter pole dominance 
fit, fixing the pole masses to the corresponding values of the vector-meson masses, m!-, 
quoted in Table 4. Diamonds correspond to the form factors at q 2 = (up to a factor Zy*) 
determined by using this second method. The fits using the constrained pole masses, also 
lead to acceptable x 2 /dof. As can be seen from these figures, both methods of extracting 
f + (0)/Zyf agree remarkably well, which gives us confidence in our procedure 11 . It can be 
also seen from these figures, that the q 2 dependence of both form factors fx(q 2 ) and f^il 2 ) 
is well described by the pole dominance model. 

For the D — > tt decay 12 we obtain, following the same steps as in the D — >■ K case: 



0.76 



.5 ; ™i 



/&(Q) 



0.74 



-4 



= 0.69 ^ ; < = 0.74 + ° [a^] 



10 This value was obtained from the analysis of scalar-scalar two-point functions, not from the study of the 
g 2 -dependence of f%{q 2 ). 

11 Because of our poor determination of the scalar + meson mass we can not do a similar comparison for 
/°(0). However, the value quoted in eq. (43) for agrees well with the value of 2.3 ± 0.2 GeV mentioned 
above ( [13]) and we expect a similar situation for /° as that obtained for / + . 

12 In order to compute q 2 and the form factors in the chiral limit we take = 0.05a -1 instead of 0, which 
corresponds to a hopping parameter of 0.14310(2) determined from eq. (24). The use of massive instead of 
massless u and d quarks, which is important for the study of the q 2 dependence of the D — > tt decay, has 
no practical consequences in the determination of the strange quark mass or in the determination of the p 
meson mass (which has been used in fixing the lattice spacing). 
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qV 



Figure 2: Results for the form factor fx(q 2 ) as a function of the dimensionless quantity 
q 2 a 2 . Left: Ki = Ki s = 0.14144. Right: Ki = K cr i t , Ki s = k s . Solid lines represent the 
pole dominance behaviour determined from a two-parameter fit to the data (parameters 
of eq. (43), for the case k\ = n cr i t , ki s = k s ). Dashed lines represent the pole dominance 
behaviour determined from a one-parameter fit to the data (fixing the pole masses to the 
corresponding values of the vector-meson masses, mi-, quoted in Table 4). Crosses and 
diamonds correspond to the values of the form factor at q 2 = (up to a factor Zy*) 
determined from two-parameter and one-parameter fits to the data, respectively. 
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Figure 3: Results for the form factor fx(q 2 ) as a function of the dimensionless quantity 
q 2 a 2 . Left: k\ = ki s = 0.14144. Right: k\ = n cr i tl Ki s = k s . The curves represent the 
pole dominance behaviour determined from a two-parameter fit to the data (parameters of 
eq. (43), for the case Ki = K cr i t , Ki 3 = k s ). Crosses correspond to the values of the form 
factor at q 2 = (up to a factor Zy f ) determined from a two-parameter fit to the data. 
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0.60 + 1 ; ; mf + = 0.91 + 13 g [a" 1 ] (44) 
/+(0) 



K 



(0) 



0.92 ± 0.08 (45) 



where we can see that the value we have obtained for the ratio of form factors /+(0) //#(0) 
is consistent with the experimental numbers quoted in section 2. Note that with the present 
errors, neither the theoretical nor the experimental result for this ratio gives clear evidence 
for iSC/(3)-fiavour violations (deviations from unity). Note also that /+(0) agrees within 
errors with /°(0) . In Figure 4 we compare the chirally extrapolated pairs (q 2 ,f^(q 2 )) with 
the pole dominance behaviour determined by the parameters of eq. (44) (solid line) and 
with the results from a one-parameter pole dominance fit, fixing the pole mass to that of 
the vector-meson mass m^L, quoted in Table 4 (dashed line). Both procedures of extracting 
/+(0) /Zy* lead again to values for /+(0) in an excellent agreement (cross and diamond in 
Figure 4). 




Figure 4: Results for the form factor f^(q 2 ) as a function of the dimensionless quantity 
q 2 a 2 . The solid line represents the pole dominance behaviour determined by the parameters 
of eq. (44). The dashed line corresponds to a one-parameter fit to the data (fixing the pole 
mass to that of the vector-meson m^L, quoted in Table 4). The cross and the diamond 
correspond to the values of the form factor at q 2 = (up to a factor Zy*) determined from 
two-parameter and one-parameter fits to the data, respectively. 
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6.2 The 0" ->■ 1" case: 



For the D — > K* decay we obtain: 



AM 



AM 

Z A 

Ao(0) 



0.67 + \ ■ mf + = 1.1 +5 [a" 1 ] (46) 
0.62 + J ; mf = 0.46 + 16 ? [a" 1 ] (47) 
0.71 + " ■ mf- = 0.59 + " [a' 1 ] (48) 



// _ 5 ) ""0- "-"^ _ 5 



V(0) 



. -1.15 : g 7 ; mf =0.85 [a" 1 ] (49) 

The results obtained for mf and m^i from the form factors V and A are consistent within 
1 or 1.5 standard deviations with the values quoted in Table 4. In the Ai case, the pole mass 
is compatible with the value of 2.5 GeV, corresponding to the D s \ resonance 13 , which was 
used in the extraction of the form factors in refs. [22]-[23]. 

In Figures 5-8 we show the form factors Ai,A 2 ,A and V as functions of q 2 for two com- 
binations of the light quark masses: Ki = Ki s = 0.14144 and = n cr it, = n s - 

In all 

cases the solid line corresponds to the comparison of the pairs (q 2 , form factor) with the 
g 2 -dependence of the form factors determined from a two-parameter pole dominance fit to 
our data (parameters of eqs. (46)-(49), for the case k\ = K criu ki s = k s ). Crosses correspond 
to the form factors at q 2 = (up to a factor Zy* or Z e ^ ) determined in this way. In 
the cases of V and A , we also compare (dashed lines) the g 2 -dependence of our data with 
that determined from a one-parameter pole dominance fit, fixing the pole masses to the 
corresponding values of the vector and pseudoscalar-meson masses, mi- and m - , quoted 
in Table 4. For the axial form factors A x and A 2 , we only make such a comparison for the 
physical situation ki = K cr i t , Ki s = k s where we fix the pole mass to the value used in [22]-[23] 
(2.5 GeV 0.9[a -1 ]). Diamonds correspond to the form factors at q 2 = (up to a factor 
Zyf or Z e /f) determined by using this second method. As can be seen from these figures, 
both methods of extracting the form factors a q 2 = agree well, serving as further check of 
consistency. 

The q 2 dependence of A and Ai is reasonably well described by the pole dominance model, 
in contrast with sum-rules calculations which predict for A\ a much weaker q 2 dependence 
than would be given by dominance of the lowest expected cs state in the J p = 1 + channel 
[37]. However, after our discussion in section 5 of the possible dependence on q 2 of the 



13 



Note however, that the spin-parity quantum numbers of this resonance have not been confirmed yet. 
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discretisation errors, we must be cautious in our affirmations about the q 2 dependence of the 
form factors and we can not draw any definitive conclusion, without a better understanding 
of the size and the q 2 dependence of the lattice artifacts present in our simulation. In the 
cases of A 2 and V, our errors are too large to determine, in a precise way, its q 2 dependence. 

As mentioned above, the values of the form factors at q 2 = have been extracted by fitting 
the chirally-extrapolated data to the pole dominance model. However, in our study we have 
some momentum channels with values of q 2 close to q 2 = 0. Another way of obtaining the 
form factors at q 2 = is to take the momentum channel which provides (in the chiral limit) 
the value of q 2 nearest to zero, and by means of the pole dominance model (with a fixed pole 
mass) extrapolate the form factor to q 2 = 0. This method for extracting the form factors 
at q 2 = has the advantage that it only requires a small extrapolation in q 2 , but on the 
negative side it only uses a single lattice point. Except in the cases of V(q 2 ) and Ai(q 2 ), the 
results obtained in such a way would agree within errors with those quoted in eqs. (43), (44) 
and (46-49). For A\ and particularly for V, the point nearest to q 2 = appears to be high 
compared to the neighbouring points, giving a higher value of the form factors at q 2 = 
if only this point is used. In the D-decay into vector mesons, the momentum channel with 
the value of q 2 nearest to zero corresponds to the transition p D = (0, 0, 0) — > \p K * \a = tt/12 
and is averaged over the six equivalent momenta of the light meson. For V(q 2 ), the other 
three channels plotted in Fig. 8, all correspond to transitions in which p D a = (1, 0, 0)7r/12. 
The difference between these two-sets of points (pu = and pr>a = (1, 0, 0)7r/12) is partly 
statistical but it may also be partly due to systematic errors affecting the two data sets 
differently. We have decided to be cautious, and to include this difference in the errors in 
our final results for Ai(Q)/Z e /^ and V(0)/Zy* . Thus, our final values for these two form 
factors are: 

0.67 + _\ (50) 
1.15 + _ 2 \ (51) 

mf + = 1.1 + _ 3 2 [cr 1 ] (52) 
mf + = 0.44 +1 [a' 1 ] (53) 
mf_ = 0.60 + _\ [a' 1 } (54) 



^(0) 



z 



eff 



no) 

z e v ff 



For the D — >■ p decay we find: 

^f(0) 



7eff 



0.60 



+ 5 
- 4 



^(0) 

ryeff 
Z A 

Z A 



0.48 



0.66 



-10 

-li 
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Figure 5: Results for the form factor Ai(q 2 ) as a function of the dimensionless quantity 
q 2 a 2 . Left: Ki = Ki s = 0.14144. Right: Ki = K cr i t , Ki s = k s . Solid lines represent the 
pole dominance behaviour determined from a two-parameter fit to the data (parameters of 
eq. (46), for the case k\ = n cr i t , ki s = k s ). The dashed line, on the right, represents the pole 
dominance behaviour determined from a one-parameter fit to the data (fixing the pole mass 
to 0.9[a _1 ] ~ 2.5 GeV). Crosses and diamond correspond to the values of the form factor at 
q 2 = (up to a factor Z e / f ) determined from two-parameter and one-parameter fits to the 
data, respectively. 
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Figure 6: Results for the form factor A 2 (g 2 ) as a function of the dimensionless quantity 
q 2 a 2 . Left: Ki = Ki s = 0.14144. Right: Ki = K cr i t , Ki s = k s . Solid lines represent the 
pole dominance behaviour determined from a two-parameter fit to the data (parameters of 
eq. (47), for the case k\ = n cr i t , ki s = k s ). The dashed line, on the right, represents the pole 
dominance behaviour determined from a one-parameter fit to the data (fixing the pole mass 
to 0.9[a _1 ] ~ 2.5 GeV). Crosses and diamond correspond to the values of the form factor at 
q 2 = (up to a factor Z e / f ) determined from two-parameter and one-parameter fits to the 
data, respectively. 
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Figure 7: Results for the form factor A (q 2 ) as a function of the dimensionless quantity 
q 2 a 2 . Left: Ki = Ki s = 0.14144. Right: Ki = K cr i t , Ki s = k s . Solid lines represent the 
pole dominance behaviour determined from a two-parameter fit to the data (parameters 
of eq. (48), for the case k\ = n cr i t , ki s = k s ). Dashed lines represent the pole dominance 
behaviour determined from a one-parameter fit to the data (fixing the pole masses to the 
corresponding values of the pseudoscalar-meson masses, m - , quoted in Table 4). Crosses 
and diamonds correspond to the values of the form factor at q 2 = (up to a factor Z e J*) 
determined from two-parameter and one-parameter fits to the data, respectively. 
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Figure 8: Results for the form factor V(q 2 ) as a function of the dimensionless quantity 
q 2 a 2 . Left: Ki = Ki s = 0.14144. Right: Ki = K cr i t , Ki s = k s . Solid lines represent the 
pole dominance behaviour determined from a two-parameter fit to the data (parameters 
of eq. (49), for the case k\ = n cr i t , ki s = k s ). Dashed lines represent the pole dominance 
behaviour determined from a one-parameter fit to the data (fixing the pole masses to the 
corresponding values of the vector-meson masses, mi-, quoted in Table 4). Crosses and 
diamonds correspond to the values of the form factor at q 2 = (up to a factor Zy*) 
determined from two-parameter and one-parameter fits to the data, respectively. 



31 



= L08 Z ; < = 0.91 « (55) 

As in the case of the D — > K* decay, we have a good determination of the A x (which domi- 
nates the decay rate) and A form factors and a poorer determination of V and A 2 . We have 
increased the upper errors of the form factors V and A\ by 0.18 and 0.01 respectively, in order 
to make our quoted values for the form factors at q 2 = compatible with the determination 
of these two form factors, from the momentum channel pr> = (0, 0, 0) — > \pk* \ = 7r/12[a -1 ]. 
In Figure 9 we compare the chirally extrapolated pairs (q 2 , A{(q 2 )) with the pole dominance 
behaviour determined by the parameters of eq. (52). 




Figure 9: Results for the form factor (q 2 , A{(q 2 )) as a function of the dimensionless quantity 
q 2 a 2 . The curve represents the pole dominance behaviour determined by the parameters 
of eq. (52). The cross corresponds to the value of the form factor at q 2 = (up to a factor 
Z e J*) determined from a two-parameter pole dominance fit to the data. 



7 Conclusions and comparison with experimental data 
and other calculations 

In this section we compare our results with the experimental measurements and other the- 
oretical calculations. In order to do this, we have to specify the values of the effective 
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renormalisation constants Z e / 3 and Z v • ; we have taken Z e / S = 1.05 _ g (eq. (42)) and 

Zy* = 0.88 ^ (eq. (41)), as discussed in section 5. In Table 10 we show our results for the 
the semi-leptonic decays D — >■ K and D — > K* . We have included the above uncertainty 
in the renormalisation constants in our final results for the form factors in order to account 
for some of the residual discretisation errors. In our quoted errors for the ratio AijA\ we 
have taken into account the fact that the discretisation errors could be different for different 
form factors and therefore in general, the effective renormalisation constant Z e J* could be 
different for A 2 than for A x giving an additional ambiguity of around 10% for this ratio (see 
eq. (42)). 

In Table 10, we also compare our predictions with the most recent experimental world average 
and with previous lattice, quark-model and sum-rule results. Our results are in reasonable 
agreement with the experimental data and the most recent lattice simulations using an 
0(a)-improved SW ([16]) and Wilson ([14]-[15]) actions. Values reported in refs. [11, 12] 
and [8, 9, 10, 13] (all of them obtained using Wilson-fermions) are also in a good agreement 
with ours, but the former are in general higher, whereas the latter are smaller, than our 
predictions. Discretisation errors are, in principle, larger for Wilson than for improved 
actions, and part of the discrepancies between different lattice results in Table 10, are due 
to different values used in the literature for the effective renormalisation constants Zy* and 

Z A ■ 

Looking now at our result for /+(0) //#(0) in eq. (45) there is no clear evidence of SU(3) 
flavour symmetry breaking and it is consistent with the experimental results. Furthermore, 
our prediction 

/+(0) =0.61 (56) 

compares well with lattice calculations obtained with Wilson fermions (0.58 ± 0.09 [10], 
0.84 ±0.12 ±0.35 [11] and 0.64 ±0.09 [15]) and other theoretical calculations (0.69 [30, 33], 
0.51 [32, 34] and 0.6-0.75 [35, 36]). The situation is similar for /°(0) . 

For the decay D — » p , we find 

yp (0) = 0.95 ; ^i(0) = 0.63 +g (57) 

^(0) = 0.51 ; ^(0) = 0.70 + J 2 (58) 

Our results are in good agreement with the most recent lattice simulation ([15]) and with the 
previous lattice calculations of refs. [10] and [12] (the results quoted in [10] are however, in 
general smaller than our predictions) and slightly smaller than the quark model prediction 
of ref. [30]. 

In the present study we have found not only the values of the form factors at q 2 = for 
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Source 






/&(0) 


V(0) 


V/A l 


Exp. 


World Ave. [28] 
World Ave. [29] 


0.77 ± 0.04 
0.70 ± 0.03 




1.16 ±0.16 


1.90 ± 0.25 


Lattice 


This work 


0.67 +7 Q 

— o 


0.65 ± 0.07 


1.01 + _Z 


1.4 + _l 


Gauge 


ELC [13] 


o.6o « +; 




0.86 ±0.24 


1.30 ±0.2 




APE [16] 


0.72 ±0.09 




1.0 ±0.2 


1.6 ±0.3 




BKS [11]- [12] 


0.90 + ; + _i 


0.70 + 8 +24 

- 8 -24 


1 43 + 45 + 48 

_ 45 -49 


1 99 + 22 +31 

^ -22 -35 




BG [15] 


0.73 ± 0.05 


0.73 ±0.04 


1.24 ±0.08 


1.79 ± 0.09 




WU [14] 


0.76 ± 0.15 


0.75 ±0.06 


1.05 ±0.33 






LMMS [8]-[10] 


0.63 ± 0.08 




0.86 ±0.10 


1.6 ±0.2 


Quark 


ISGW [32] 


0.76 - 0.82 




1.1 


1.4 ± 0.4 


Models 


WSB [30] 


0.76 




1.27 


1.4 




KS [33] 


0.76 




0.8 


1.0 




GS [34] 


0.69 




1.5 


2.0 


Sum Rules 


BBD [37] 

AEK [35] 
DP [36] 


0.60 +Jjj 

0.60 ± 0.15 
0.75 ±0.05 




1.10 ±0.25 


2.2 ±0.2 


Source 




^(0) 


A 2 (0) 


A 2 /A x 


A 


Exp. 


World Ave. [28] 


0.61 ± 0.05 


0.45 ±0.09 


0.74 ±0.15 




Lattice 


This work 


0.70 ^ 


o.66 + ;° 

— 15 


0.9 ±0.2 


0.75 +^ 


Gauge 


ELC [13] 


0.64 ± 0.16 


0.40 ±0.28 ±0.04 


0.6 ±0.3 






APE [16] 


0.64 ± 0.11 


0.46 ±0.34 


0.7 ±0.4 






BKS [12] 


n + 14 + 28 

- 83 -14 -28 


U - Dy -14 -23 


0.70 +16 +20 

-16 -15 


94 +9 +22 

-9 -24 




BG [15] 


0.66 ± 0.03 


0.42 ±0.17 


0.71 ±0.20 






WU [14] 


0.59 ± 0.08 


0.56 ±0.40 








LMMS [8]-[10] 


0.53 ± 0.03 


0.19 ±0.21 


0.4 ±0.4 




Quark 


ISGW [32] 


0.8 


0.8 


1.0 ±0.3 




Models 


WSB [30] 


0.88 


1.15 


1.3 






KS [33] 


0.82 


0.8 


1.0 






GS [34] 


0.73 


0.55 


0.8 




Sum Rules 


BBD [37] 


0.50 ± 0.15 


0.60 ±0.15 


1.2 ±0.2 





Table 10: Form factors at q 2 = for the semi-leptonic decays D — >■ K and D — > K* : 

comparison of our results with experimental data and with other theoretical calculations. In 

obtaining our results we have used Z e /* = 1.05 ^ and Zy* = 0.88 All lattice gauge 

calculations have been obtained using Wilson fermions except that of ref. [16] and the present 

work, where an 0(a)-improved SW-action has been used. 
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the different decay processes, but also their q 2 dependence in a wide region around q 2 = 0. 
Thus we can estimate the integrals of eqs. (8-9) and obtain the total decay rates. Our 
predictions, together with the experimental measurements and other theoretical calculations, 
are presented in Table 11. As mentioned above, we have a poor determination of the q 2 
dependence of the A 2 and V form factors in the D decays into vector mesons. However, 
as can be seen in eq. (9), the contribution of these form factors to the decay rates is small 
and only important in the proximity of q 2 = 0, where their contributions are reasonably well 
determined. For example, pole masses, for the A 2 form factor in the D — > K*, p decays, three 
times larger than those of eqs. (47,53) give total decay rates and ratios T L /T T which differ 
from those quoted in Table 11 only at the level of (0.3-0.5) standard deviations. Therefore, 
we are confident that we can use the values of the form factors from our simulation for 
calculating the total decay rates. 

In the 0~ — > 0~ case, (Zy^) 2 is an overall factor in the expression for the width and thus 
we could quote our result for the decay rate in terms of (Z^) 2 /0.88 2 . However, in the 
decay into vector mesons, the form factors H ± mix the contribution of both the vector and 
the axial form factors and thus such a factorization cannot be made. Therefore in both 
cases (decays into pseudoscalar and vector mesons) we have decided to include in the quoted 
statistical errors of our results, the uncertainty due to Zy* and Z e J* . We have estimated 
this uncertainty by computing the extreme values which would be obtained for the different 
decay rates if the errors of eqs. (41-42) were taken into account. On the other hand, the 
q 2 dependence of the form factors is determined by the different pole masses, quoted in 
eqs. (43-55), whose physical values depend on the precise value taken for the lattice spacing, 
a -1 , and thus the results obtained for the decay rates will also depend on the scale a -1 . The 
second set of errors in our results of Table 11 is due to the uncertainty in the determination 
of the lattice spacing; we have taken a -1 = 2.85 ± 0.15 GeV. This ambiguity in the scale 
has, in general, a small effect on the decay rates, and in some cases is negligible. 

As can be seen in Table 11, our results are in excellent agreement with the experimental 
data. This agreement, together with that already shown in Table 10, provides further con- 
fidence that lattice QCD is becoming a reliable quantitative tool for non-perturbative QCD 
phenomenology. Studies of charm physics on the larger lattices which will shortly become 
available, will provide a fruitful area of investigation, and will enable the control of the sys- 
tematic errors (except quenching) present in these calculations. This understanding of the 
discretisation errors will make it possible to obtain accurate estimates of the the QCD-non 
perturbative corrections to the B — >■ n and B — ?► p decays, from which we expect to extract 
the Cabibbo-Kobayashi-Maskawa matrix element \V u b\. 

We end this paper with a brief summary. The study presented in this paper is one of the first 
calculations of the form factors of weak vector and axial currents relevant for semi-leptonic 
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Source 




T(D ->• K) 


r(D ->• k*) 




T(D p±) 


Exp. 


World Ave. [28] 

Wnrlrl Aire T9Ql 

vvoiici rvve.[^yj 


9.0 ±0.5 

1 . 1 Hz U.D 


5.1 ±0.5 


0.60 ±0.15 




Lattice 
Gauge 


This work 

ELC [13] 
APE [16] 
LMMS [10] 


7.0 ± 1.6 ±0.4 

5.4 ±3.0 ±1.4 
7.8 ±2.2 
5.8 ± 1.5 


6 - -1.6 

6.4 ± 2.8 
6.3 ± 1.7 
5.0 ±0.9 


0.52 ±0.18 ±0.04 

0.5 ±0.3 ±0.1 
0.7 ±0.2 
0.5 ±0.2 


0.43 ± 0.11 

0.60 ±0.3 ±0.1 
0.5 ±0.2 
0.40 ± 0.09 


Quark 
Models 


ISGW [32] 
WSB [30] 
KS [33] 
GS [34] 


8.5 
8.26 
10.2(e) - 9.9(/i) 
7.1 


9.1 ±0.25 






Sum 
Rules 


BBD [37] 
AEK [35] 
DP [36] 


6.4 ± 1.4 

5.1 ± 1.7 

8.2 ± 1.1 


3.2 ± 1.3 


0.76 ±0.24 





Source 




V{D^K*) 


(ft)- 




Exp. 


World Ave. [29] 


0.55 ±0.07 


1.2 ±0.1 




Lattice 


This work 


0.86 ±0.28 ±0.03 


1.06 ±0.16 ±0.02 


1.05 !^±0.04 


Gauge 


ELC [13] 


1.1 ±0.6 ±0.3 


1.4 ±0.3 






APE [16] 


0.8 ±0.3 


1.3 ±0.3 






LMMS [10] 


0.86 ±0.22 


1.51 ±0.27 


1.86 ±0.56 


Quark 


ISGW [32] 


1.1 


1.1 ±0.2 




Models 


WSB [30] 


1.15 


0.9 






KS [33] 


0.95 


1.1 






GS [34] 


1.4 


1.2 




Sum Rules 


BBD [37] 


0.5 ±0.15 


0.86 ±0.06 





Table 11: Semi-leptonic partial widths for D — > K,K*,tt and p, using \V CS \ = 0.975 and 
\V c d\ = 0.222. We also report the ratio of the longitudinal to transverse polarisation partial 
widths for D ->■ K* and D -)> p. Units in 10 10 s _1 . Ref. [28] gives T(D tt) = 1.2 ±0.3. We 
have assumed isospin symmetry and we have taken a value for the decay rate, with charged 
pions in the final state, of T(D — > 7r ± ) = 0.60 ± 0.15. 
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decays of D-mesons, performed using the improved quark action proposed by Sheikholeslami 
and Wohlert [7]. Our results for the form factors (Table 10) and decay rates (Table 11) have 
reasonably small errors and are in good agreement with experimental measurements. The 
results at non-zero momentum transfer are, in general, in agreement with the pole dominance 
model. 

We have tried to minimize systematics by working with an improved action to reduce dis- 
cretization errors, and on fairly large volume in the hope that finite-size effects would be 
small. Nevertheless, it is important that our simulation be repeated on lattices of different 
sizes and spacings in order to quantify more precisely the systematic effects, which could 
modify the results presented in this work, in particular the g 2 -dependence found for the 
different form factors. 
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